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A new transform method for solving boundary value problems for linear and integrable nonlinear 
PDEs recently introduced in the literature is used here to obtain the solution of the modified 
Helmholtz equation qxx{x, y) + qyy{x, y) — Aj3^q{x, j/) = in the triangular domain Q<x<L~y< 
L, with mixed boundary conditions. This solution is applied to the problem of diffusion-limited 
coalescence, A -\- A ^ A, in the segment (— L/2, L/2), with traps at the edges. 
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I. INTRODUCTION 

A new method for solving boundary value problems for linear and for integrable nonlinear PDEs has been introduced 
recently |^ . Here we apply this method to the equation 

E^^+Eyy+j{~E,+Ey)=Q, (l.l) 

in the triangular domain —L/2 < x < y < L/2, where E{x,y) is a scalar function and 7 is a positive constant. 
A solution of Eq. (1.1) in the semi- infinite wedge < x < ?/ has been presented in Using the substitution 



E{x,y) = 1 — e 2 ~^q[x,y), Eq. (1.1) becomes the modified Helmholtz equation 



q.. + qyy - = , /5 = ^ • (1-2) 

Eq. ( pTll ) with ^ — v/D represents the steady state of the diffusion- limited reaction A-\- Aotl the line, where the 
A-particles diffuse with diffusion constant D, they merge immediately upon encounter, and split into two particles (the 
back reaction) at rate v |^-|^. E{x,y) represents the probability that the interval {x,y) is empty. The concentration 
profile of the particles is related to E[x, y) thr oug h c{x) ~ —Ey{x, x). Suppose that we limit ourselves to the segment 



—L/2 < x < L/2, then the domain of Eq. (Id) is —L/2 < x < y < L/2. The forward reaction is described by 



the boundary condition (BC) E{x,x) = 1. If there are perfect traps at the edges, x — ±L/2, one gets the BCs 



E^{-L/2,y) = 0, Ey{x,L/2) = 0. These BCs transform into the following BCs for Eq. (|lj) 

q{x,x)=0, -^<x<^, 

i'7(-iy) + fe(-iy) = 0, -|<y<-|, 

-^qix,^) + q.,ix,^)^0 , -^<x<^. 

We rotate and translate the (x,j/)-axes, with the mapping {x,y) t-^ {—y + L/2,x-\-L/2). Eq. (|l.2| ) remains invariant, 
but the domain is now 0<x<L — y<L — the isosceles right triangle with vertices at (0, 0), (0, L), {L, 0) — and 
the BCs become 

q{x,L-x) = , 0<x<L, (1.3a) 
^q{x,0)-^qy{x,0) = , 0<x<L, (1.3b) 
^q{0,y) + qAO,y) = 0, 0<y<L. (1.3c) 



For the sake of generality, instead of the BC ( |l.3| b), (1.3c) we consider 
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qy{x,0) 



fiy), 



0<x<L, 

o<y<L, 



(1.3b') 



(1.3c' 



where /(•) is an arbitrary smooth function. 

We win show that: (a) Eg. (1.2) with the BCs ( |1.3| a), (|l.3| b'), (P^ c') has a unique solution that can be expressed in 
closed form, (b) Eq. (O) with the homogeneous BCs (|1.3|a), (|1.3| b), (|l.3|c) has only the trivial solution q{x,y) = 0, 
i.e., the only steady state of the process A-\- A ^ A, in a segment demarcated by traps, is the vacuum — when there 
are no particles left — regardless of the magnitude of v, the rate of the back reaction A A + A. (c) For large 
back reaction rates, 7L 3> 1, the characteristic relaxation time to the empty, absorbing state grows exponentially as 

2\„vL/2D 

+ iy, let a bar denote the complex conjugate (z = 
y < L (see Figure 1); 



{D/2v)e 
Let z = 
0< X < L 



iy), and let Zj denote the corners of the domain 



Z2 



0, 



Z3 



(1.4) 



II. THE ANALYSIS OF THE INHOMOGENEOUS PROBLEM 

It is shown in |^ that the general solution of the modified Helmholtz equation in the above domain can be represented 



as 



q{x,y) 



2m 



E 



0<x<L-y <L 



(2.1) 



where £1, £2, £3, are the rays on the complex fc-plane defined by argfc = 0, ^, and oriented from zero to infinity 
(see Figure 2), while the functions Pj{k) are defined by 



p,{k) 



-ikz+il-r-zf 1 



(qx - iqy)dz + i^qdz] 



fc e c 



1,2,3, 



Z4 — Zl 



Using the boundary conditions (1.3) to simplify the expressions for Pj{k), we find the following: 

pi(fc) = -ig(0,0) + ia(fc)V'i(-ifc)-«F(-?;fc) , keC , 
P2{k) = ^9(0, 0) + ia(-ifc)V'2(fc) - iF{k) , fc £ C , 

P3(fc) =i^(fc)V'3(-fce^^) , fceC, 

where 



a(fc) = i(^ + fc + i), 



E{k) 



2 7o 



e(''+^)yf{y)dy. 



(2.2) 

(2.3a) 
(2.3b) 
(2.3c) 

(2.4) 



and the unknown functions "01, ip2, tp3 are defined by 



Mk) 



^''+^^^q{x,0)dx , Mk)^ f e(''+^'>yq{0,y)dy , Mk) ^ 



V2L 



qsi 



V2' 



V2 



ds 



Indeed, for the derivation of ( |2.3| a) we use z = x, and we note that the boundary condition ( |1.3| 3') implies 
i(g,(x, 0) - iqy{x, 0)) + i^qix, 0) = ig,(x, 0) + + ^)q{x, 0) - |/(x) ; 



(2.5) 



integrating by parts the terms involving 173; we find ( |2.3| a). The derivation of ( |2.3| b) is similar, where we use the 
condition ( |l.3| c'). For the derivation of ( p.3p ) we use z = iL + x — ix, and we note that the boundary condition 
q{x, L — x) — implies q^ix, L — x) — qy{x, L ~ x) — 0. 

In order to simplify the analysis, we have assumed that the same function / appears in the BCs ( [l.3| 3') and ( |l.3| c'). 
This implies that the PDE (1.2), the triangular domain, and the BCs (1.3a), ([l.Spj'), (pT^c') are invariant under the 
reflection a; <-> y, thus q{x,y) — q{y,x). Hence, 4'i{k) — ip2{k). 

We introduce the following notations: 



Mk)=Mk)=m, M-ke-'^'^)=^{^) , i^3(-fce-/4) = 0(+) , e(fc, z, z) = e'^^"^"^^ 



(2.6) 
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A. The Analysis of the Global Relation 



Eqs. (2.3) express Pj{k) in terms of the unknown functions ip{—ik), (p{k), and These functions satisfy the 



global condition: X]j=i Pji^) — lO- This equation, and its complex conjugate, are 

aik)ip{-ik) + a{-ik)ip{k) + E{k)iP{+) = F{k) + F{-ik) , keC 
a{k)Lp{ik) + a{ik)ip{k) + E{k)i:{-) = F{k) + F{ik) , keC . 



(2.7) 
(2.8) 



Following ||6| we supplement these equations with the equations obtained from Eqs. ( |2.7[ ), ( |2.8| ) by using the 
transformations in the complex fc-plane which leave invariant the pairs {ip{—ik), ip(ik)} and The first 

pair is invariant unde r fc i—^ — fc, a nd the second pair is invariant under {k t-^ —ik,k t—f ik}. Using the latter 
transformations, Eqs. (^) and (^) yield 



a{-ik)(p{~k) + a{-k)ip{-ik) + E{-ik)il;{-) = F{-ik) + F{-k) , keC 
a{ik)ip{-k) + a{~k)Lp{ik) + E{ik)^P{+) = F{ik) + F{-k) , fc e C . 



(2.9) 
(2.10) 



Eqs. (2.7)- (|2.10| ) are invariant under fc i— > — fc, thus we do not obtain any additional equations using this transforma- 
tion. Eqs. ( |2.7] )-(2.10) are the basic equations needed for the determination of the unknown functions ^pik), (^(— ifc), 
'(/'(-I-). The analysis of the basic equations leads to a matrix Riemann-Hilbert problem. However, in what follows we 
will show that this problem can be bypassed, and that q{x, y) can be obtained using only algebraic manipulations of 
the basi c eq u ation s. 

Eqs. (2.7)-( 2.1C ) imply that Lp{~ik), (p{k), can be expressed in terms of Lp{ik) and 



ip{~ik) = A{—ik)ip{ik) + 



a(-fc)A(fc) 



[A{-ikfE{ik) ~ E{~ik)]^|;{-) + Gi(fc) 



m = -^M^k) - ^M-) + G2(fc) , 



a{ik) 



a{ik) 



where 



^ A{-ik)E{k)^A{k)E{~ik) 
^ ' A(fc) 

a(k) 



G3(fc) 



m 



a{~k) 



A{k) = E{k) + A{k)A{~ik)E{ik) 



and the known functions Gj{k), j — 1, 2, 3, are defined in terms of / as follows: 



G'i(fc) 



A(fc)Q!(-fc) 



{[E{k) + A{-ik)E{ik)][F{-ik) - A{~ik)F{ik)] 
+ [1 - A{-ik)][A{-ik)E{ik)F{k) + E{k)F{-k)]} . 



(2.11a) 
(2.11b) 
(2.11c) 

(2.12) 
(2.13a) 



G2(fc) 



F{k)+F{ik) 
a{ik) 



(2.13b) 



G3(fc) 



A(fc) 



{[1 - A{k)][F{-ik) - A{-ik)F{ik)] + [1 - A{-ik)][F{k) - A{k)F{-k)]} 



Indeed, Eq. ( |2.1l| b) is Eq. (gj). Eliminating (p{-k) from Eqs. (|2J), ( p.lOD , we find 



a{-k)ip{-ik)+E{~ik)ilj{-) - A{-ik)[a{-k)ip{ik) + E{ik)ilj{+)] 
^ F{-ik) + F{-k) - A{~ik)[F{ik) + F{-k)] . 



(2.13c) 



(2.14) 



Replacing in this equ ation (p(ik) by Eq . (|2.8D and comparin g wit h Eq. (2.7) we find Eq. ( 2.11| c). Replacing V'(+) in 
terms of ^{-) in Eq. ( |2.14| ), using Eq. ( |2Tl| c), we find Eq. ( ^.ll|a,). 

Eq. ( ^.l[ ) expresses q{x, y) in terms of Pj{k), and Eqs. ( [2.3| ), (2.11) express Pj{k) in terms of the unknown functions 
fiik), ip{—), and the known functions Gj{k). The known functions give rise to the contribution 



3 



G{x,y) - i-^ e{k,z,z)[a{k)G,{k) - F{-ik)\^ + 



dk 

e(fc, z,z)[a{~ik)G2{k) - F{k)] — 

k 



1 f dk 

-—J^eik,z,z)E{k)G^{k)^. (2.15) 



In what follows wc will show that, by using appropriate contour rotations, the integrals involving the functions ip{ik), 
ip{—) can be evaluated in terms of residues. Furthermore, these residues can be computed in terms of the functions 
Gj{k). For the justification of these rotations we use the following facts (see Figure 3). 

• e(fc,z,j), e{k, z,z)E{k), e{k, z,z)E{~ik), are bounded for < argfc < f j f < argfc < < arg/c < 27r, 
respectively. 

• E{-k)E{ik) and ip{-) are bounded for < arg/c < while E{k)E{-ik)ip{-) is bounded for ^ < arg/c < ^. 

• A(/c) ~ E{k), /c ^ and /c oo, in -f < argfc < ^; A(fc) - E{ik), fc and fc ^ cx), in ^ < argfc < ^. 

Indeed, since a; > and y > 0, e(fc,z,z) is bounded both at fc = and fc = oo in the first quadrant of the 
complex fc-plane. Since — f < arg(z — 23) < — f , it follows that if ^ < argfc < ^ then < argfc(z — Z3) < tt. 

Hence exp[ifc(2; — 2:3) — i^{z — 03)] is bounded both at fc = and fc = oo; using Z3 = this exponential equals 
e{k, z,z)E{k). Similar considerations apply to e{k, z,z)E{—ik). 

A(fc) = E{ik)[E{k)E{~ik) + A{k)A{-ik)]. If -f < argfc < ^, E{k)E{-ik) is exponentially large at fc = 
and fc = CX3, and A(fc) ^ i?(fc). Similarly, if < argfc < E{k)E{—ik) is exponentially small, and A(fc) ~ 
£:(ifc)^(fc)A(-ifc) - E{ik). 

involves -(fce"*^/'' + ^e"/"), thus it is bounded for -f < argfc < Similarly for E{k)E{-ik)^:{-). 
The contribution of the integral along £3, due to the terms involving (see Eq. ( ^.ll| c)), gives rise to two 

integrals: one involving e{k,z,^)A{k)E{—ik)E{k)'4'{—)/kA{k), and one involving e{k, z,^)A{—ik)E{k)'^ /kA{k). The 
first integral is bounded in ^ < argfc < 2tt, while the second integral is bounded in | < argfc < Indeed, the 
integrand of the first integral is dominated by 

[eik,z,z)E{-ik)][E{k)E{~-ik)ij{-)] , 5il<argfc<^; [e(fc, z, z)^(-ifc)][^(-)] , ^ < argfc < 2^ , 

and each of the brackets is bounded. Similarly, the integrand of the second integral is dominated by 

[eik,z,z)E{k)][ip{-)] , |<argfc<^; [eik,z,z)E{k)][E{k)E{-ik)^P{-)], f < argfc < ^ , 

and each of the brackets is bounded. 

Hence, the integral along €3, due to the terms involving ipi—), equals an integral along £1 involving 
e{k,z,T)A{k)E{—ik)E{k)ip{—)/kA{k), an integral along £2 in volvin g e (fc, z, j)A{—ik)E{k)'^ /kA{k), and a contribu- 
tion due to residues which will be computed below (see E qs. (2T9| b), ( 2.19| c)). Combining these integrals with the 
integrals due to (p{—ik) and to (p{k) (see Eqs. (2.11a) and (2.11b)), we find 



1 /■ dk 

Ji{x,y)^— j e{k,z,z)[^q{Q,Q) + a{k)A{-ik)^{ik)]— , (2.16) 

Ux,y)^^ j ^^^^A{k)A{-^krE{k)E{^k)^P{-)f. (2.17) 

For fc in the first quadrant of the complex fc-plane, E{k) / A(k) is dominated by 1, and each of the terms e{k,z,'z), 
ip{ik), E{ik), V-'(— ) is bounded. Thus, both Ji and J2 can be computed in terms of residues. 
The definition of ^(fc) implies 

■,^ (fc + Ai)(fc + A2) ,,,, (fc-FAi)(fc- A2) , 7,, , 7, 

-^(--^)-- (fc-A0(fc-A2r ^(^^-- (fc-A0(fc+A2r ^^ = i(-i+^)' 

so the poles of A{—ik) and A{k) occur at Ai, A2, and at Ai, — A2, respectively. Similarly, the poles of A{ik) and 
A{—k) occur at — Ai, — A2, and — Ai, A2, respectively. Using these facts it follows that 

3 

q{x,y) = G{x,y) +J2Rjix,y) + Pi^.y) , (2-18) 



4 



where G{x,y) is defined by Eq. (2.15), P{x,y) is the contribution to Ji and J2 due to the poles of a{k)A{—ik) and 
of A{k)A{—ik)'^, and the Rj are defined as foUows: 



3 



R2=-^Y: 



e{Kj, z, z)A{Kj)A{—iKj)'^E{Kj)E{iKj) 
KjA'(kj) 

e{X„z,z)A{-X,)E{X,f 



A,A'(A,) 



A(A2) 



R3 



e(-A2, z, z) 
A(-A2) 



E{-A2)E{iA2M-^2) 



ipik), A'(fc) denote 



^l;{k) = V(-fce-*''/'') , A'(fc) = 



dA{k) 
dk 



(2.19a) 
(2.19b) 
(2.19c) 

(2.20) 



and KjjXj, IX j denote the zeros of A(fc) in < argfc < f , f < argfc < ^ < arg/c < 27r, respectively. Multiplying 
Eq. (2.II3) by A(fc) and evaluating the resulting expression at kj — {kj, Aj, /ij}, we find 

_ [1 - A{k,)][F{-ik,) - A{-ik,)F{zk,)] + [1 - A(-»fc,)][F(fc,-) - yl(fc,-)F(-fc,)] ^ ^(^^.^ _^ q ^ (2.21) 



where 



<5(fc,) 



5(k) = A{-ik)E{k) + A{k)E{-ik) , 



(2.22) 



Noting that a(fc) = (fc + Ai)(fc - A2]/2fc, a(ifc) = -(fc - Ai)(fc - A2)/2ifc, and evaluating Eq. ( |2.8| ) at fc = A2, we find 
'0(A2). Similarly, evaluating Eq. (2.9) at fc = — A2 wc find ■(/;(— A2): 



V'(A2) = £;(-A2)[F(A2) + F(zA2)] , ^^(-Aa) = E{-iK2)[F{k2) + ^(^2)] 



(2.23) 



The term P{x,y) arises from a{k)A{—ik) in Ji, and A{k)A{—ik)'^/A{k) [A(fc) 7^ 0] in J2, each of which has a 
simple pole at A: = Ai . Evaluation of the pertaining residues yields 

P{x,y) = 2e(Ai,z,z)[a(Ai)(/.(iAi) + £;(Ai)^(Ai)] . 

Evaluating Eq. ( |2.8| ) at fc = Ai, we find 

a(Al)(^(^Al) + £;(Ai)^(Ai) - F{Ai) + F(iAi) . 

Thus, 

P{x, y) = 2e(Ai, z,z)[F{Ki) + F(iAi)] . (2.24) 



In summarv : Assume that 5{kj) ^ 0, wh ere kj is a zero 0/ A(fc), and S{k), A( fc) ar e defined by Eqs. (2. lib), (2.2i), 
respectively. Then q{x, y) is given by Eq. ( 2. IS ), wh ere G{x, y) is defined by Eq. P{x, y) is defined b y Eq. ( 2.24 ): 

and Rj{x,y), j = 1,2,3, are defined by Eqs. (2.1!\), with ij{kj), '(/'(A2), '(/'(-A2) defined by Eqs. (Cil), 



III. THE PHYSICAL PROBLEM 



The physical problem corresponds to the homogeneous BCs, i.e., / = 0. In this case Eq. ( 2.1^ ) yields q{x,y) — 



Thus, we only need to consider the assumption that 5{k.n) 7^ 0. If this assumption is violated then the equations 
A(fc„) = and 5(fc„) = can be rewritten in the form 

A{iknfE{~ikr,f = 1 , (3.1a) 
A{-k^fE{k^f = 1 . (3.1b) 

Eqs. ( [3.1[ ) do not have a solution f or g eneric values of 7. Indeed, consider first the limit of infinite back reaction 
rate, 7L 00. Inspection of Eqs. ( ^ ) in this limit yields the asymptotic solution fcoo = ±Ai, ±A2. If there exists 



5 



a steady state other than q{x, y) = 0, then it would also exist for 7L large but fini te. We therefore seek solutions 
of (^) of the form fc = fcoo + e. Such solutions do not exist: Using koo = Ai, Eq. (3Ja) yields e = ^(—1 + «) 



to first order in e — while Eq. dS.lt p) yields the contradictory result e = -^{+^ — *)• The other values of koo lead to 
similar contradictory results. Thus, the only solution to the physical problem is q{x^ y) — 0, which corresponds to the 
trivial case of the vacuum; when no particles are left in the system. 



Finally, consider the relaxation of the system into the absorbing empty state. Instead of Eq. (1.1), we need to study 

E'xx + Eyy + ^{—Ex + Ey) — Et , (3-2) 
where r — Dt is a rescaled time parameter. We turn this into an eigenvalue problem, by writing E{x,y,t) = 



1 — e °"^e aCi* '■^^q„{x,y). This results in an equation for identical to Eq. (L2), valid over the same domain, but 



with 4/3 = ^7 — a. The BCs for this equation are identical to (|l.3|). We have already seen that the problem admits 
no zero eigenvalue: qo{x, y) = 0. The analysis for cr > proceeds along the same lines. Once again, the critical issue is 
whether there exist solutions of Eqs. (3.1). This time the asymptotic solution for 7L — > 00 is koo — 4 [ili*(l~4^)^^^]i 



j[±(l — 4^)^/^±i]. A perturbation analysis shows that solutions exist for finite 7L ^ 1, provided that a ~ 27^6 ^^Z^. 

The relaxation time to the empty state is therefore {Da)~^ — {D /2v^)e^'"/'^^ . 

It is instructive to compare our analysis of A + A ^ A to the mean- field result. The reaction-diffusion equation for 
the steady state of the process, in a segment demarcated by traps, is 

Dpxx + kip-k2p^ , -L/2<x<L/2, (3.3) 

where p{x) is the local particle density, fci is the rate of the back reaction A ^ A + A, k2 is the rate oi A + A ^ A, 
and the traps impose the BCs p{±L/2) — 0. This equation predicts a transition from an empty state (p — 0) to an 
active state {p > 0), when ki exceeds a certain critical value 0,^. Our exact analysis shows that in the actual system 
of one-dimensional coalescence the noise destroys the transition and the only existing steady state is the empty state. 
The nontrivial steady state of the mean-field case is echoed in the exponentially large relaxation time found for large 
back reaction rates. Although the lack of a transition cannot be established from numerical simulations, especially in 
view of the long relaxation times for jL large, previous work had suggested that a transition does not take place . 
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FIG. 1. Domain of the modified Helmlioltz equation, Eq. ( [l.2| ) 
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FIG. 2. The rays Ij, in the complex plane, along which q{x,y) is computed (Eq. (2.1)) 
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e(k,z,z)E(k) 



e(k,z,z) 




e(k,z,z)E(-ik) 



A(k)^E(ik) 
E(k)E(-ik)^(-) 



A(k)~E(k) 
E(-k)E(ik) 



FIG. 3. Regions where e{k,z,z), e{k, z, z)E(k), e{k, z, z)E{—ik), i^^—), E{k)E{—ik)^{~) are bounded, and dominant be- 
havior of A(fc). 
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